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Abstract. The evolution of the = 28 shell closure is investigated far from stability. 
Using the latest results obtained from various experimental techniques, we discuss the 
main properties of the A^ = 28 isotones, as well as those of the N = 27 and A^ = 29 
isotones. Experimental results are confronted to various theoretical predictions. These 
studies pinpoint the effects of several terms of the nucleon-nucleon interaction, such as 
the central, the spin-orbit, the tensor and the three-body force components, to account 
for the modification of the A^ = 28 shell gap and spin-orbit splittings. Analogies 
between the evolution of the N = 28 shell closure and other magic numbers originating 
from the spin-orbit interaction are proposed (A^ = 14,50,82 and 90). More generally, 
questions related to the evolution of nuclear forces towards the drip-line, in bubble 
nuclei, and for nuclei involved in the r-process nucleosynthesis are proposed and 
discussed. 

PACS numbers: 21.10.-k, 21.30.-x, 21.60.-n, 27.30.+t, 27.40.-l-z 
1. Introduction 

Guided by the existence of a spin-orbit coupling in atomic physics, M. Goeppert 
Mayer [1] and O. Haxel et al [2J independently postulated in 1949 the existence of 
a strong spin-orbit (SO) force in the one-body nuclear potential to reproduce the so- 
called 'magic numbers' above 20. Their SO interaction has to be attractive for nucleons 
having their angular momentum aligned with respect to their spin (denoted as i^) and 
repulsive in case of anti-alignment (£4). Its strength was constrained to reproduce the 
size of the shell gaps at 28, 50, 82, and 126 derived from experimental data obtained in 
the early days of nuclear structure studies. 

During the last two decades, major experimental and theoretical breakthroughs 
have strengthened that, though this simplified mean field picture of the SO interaction 
is elegant and simple, the observed SO splitting comes from the combination of complex 
effects due to the nuclear force. In particular, the size of the = 28 shell gap is 
governed by the spin-orbit, the tensor, the three-body force components, and possibly 
the change of nuclear forces at the drip-line due to interactions with the continuum. 
For instance, the A^ = 28 shell gap grows by about 2.7 MeV between ^^Ca and ^^Ca to 
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reach a value of 4.8 MeV. This effect has recently been mainly inferred to three-body 
forces. Then, starting from the doubly-magic lo^ajg nucleus and removing protons, 
a progressive onset of deformation is found to occur: igAr is a moderately collective 
vibrator, fgS exhibits a prolate-spherical shape coexistence, the ^481 nucleus is a well 
oblate nucleus, and the newly discovered f^^g is likely to be a prolate nucleus. Such 
a progressive change of structure is triggered by the combination of three effects: the 
reduction of the spherical = 28 shell gap, the probable reduction of the Z = 14 
sub-shell gap and the increase of quadrupole excitations across these gap, i.e. on the 
one side between the neutron I/7/2 and 2^3/2 orbits and on the other side between the 
proton l(i5/2 and 2si/2 orbits which are separated by two units of angular momenta {2h). 

This smooth but continuous change of nuclear shapes at the semi-magic = 28 
nuclei is quite unique in the chart of nuclides. Therefore from this study we can learn 
how a spherical rigid system could evolve towards a very collective system when a few 
nucleons are removed. There, proton-neutron forces act between nucleons in orbits 
which are different from the ones occupied in the valley of stability, hereby involving 
new facets of the nuclear forces which were not yet tested. 

The magic number 28 is in fact the third number of nucleons which is created by the 
SO splitting the first two SO gaps being 6 and 14. Similar components of the nuclear 
forces are expected to play a role in the gap values of all the SO magic numbers, leading 
to possibly similar consequences throughout the chart of nuclides. The magic number 28 
is a target choice to study nuclear forces as it involves nuclei having an intermediate mass 
and size, with orbits relatively well separated from the neighboring ones. This feature 
is essential to distinguish the effects of nuclear forces between nucleons occupying orbits 
having different properties (such as the number of nodes, the angular momentum and 
the spin orientation). Conversely as the size of nuclei grows, the in- medium nuclear 
force is weaker and modification of shell structures occur only when many nucleons 
are involved. Added to this, the nuclear orbits are more compressed in energy and 
their occupancy is diluted by nucleon-pair scatterings. It follows that the action of 
specific nuclear forces between nucleons in well established orbits is much more difficult 
to disentangle in heavy than in light systems. At the other extreme, if the system is 
composed of a small number of nucleons, the change of only few nucleons leads to very 
rapid and drastic structural changes and the very concept of the mean-field approach 
becomes inappropriate. To quote one example of quick change of a SO magic gap, the 
A^ = 14 gap collapses already after the removal of two protons, from ^^O to ^°C. 

On the experimental point of view, an incredible wealth of information has been 
obtained these last two decades by studying the evolution of nuclear structure south to 
the doubly-magic 2oCa28 nucleus. Different experimental techniques are used worldwide, 
such as P or isomeric-decay studies, atomic mass measurements. Coulomb excitation, 
transfer reaction, or in-beam 7-ray spectroscopy. These studies became feasible since 

I Note that from now on, we will rather explicitly mention SO splitting rather than SO interaction 
as the shell gaps are created by several components of the nuclear force, and not exclusively by the 
spin-orbit force. 
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Figure 1. Part of the chart of the nuchdes which wih be considered in the present 
work. Stable nuclei are displayed with a dark blue color. In the bottom part, the 
orbits of the proton and neutron valence spaces under study are shown. 



increasing beam intensities of '^^Ca and detectors capabilities become available. With 
8 neutrons more than the most abundant ^°Ca isotope, the ^®Ca isotope represents a 
fraction about 0.2% of the total Ca isotopes on earth. However it is more than 50 times 
more abundant than the (A^ — 2) nucleus ^^Ca on earth, a feature which is not so common 
in the chart of nuclides. This is ascribed to its double magicity which makes it produced 
in significant amount and poorly destroyed in specific stellar explosive environments. In 
some sense, nature has offered with the neutron-rich and doubly magic ^^Ca isotope a 
rather unique tool to study the behavior of nuclear forces far from stability. 

We shall explore in the next Sections the nuclear structure evolution of the nuclei 
presented in Figure [T] In Section |2} the general properties of atomic masses, first 
excited states and reduced transition probabilities will be presented. In Section [3} we 
shall discuss upon the evolution of proton orbits when moving from = 20 to = 28 
by studying the odd-Z isotopic chains of 19K, 17CI and 15P. Sections |4] and |5] show the 
evolution of the = 28 gap from the behaviors of the N = 27 and A^ = 29 isotones. 
Some experimental and theoretical challenges related to the nuclear forces responsible 
for the SO magic numbers will be proposed in the last Section. 

Sections 2 to 5 of the present paper update our previous work on the evolution of 
the A^ = 28 shell closure [3]. 
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Figure 2. Evolution of tlie two-neutron separation energy in the even-Z isotopic 
chains (experimental atomic masses from Refs. [HIS]). The sharp drop in observed 
in the Ca isotopic chain after = 28 progressively disappears far from stability. 

2. View of the shell-structure evolution from atomic masses, E(2f ,4:^) 

and B{E2; 0^" -> 2^) values 

2.1. Atomic masses 

The two-neutron separation energy S2n{A) can be obtained in a given isotopic chain 
from the atomic-mass differences between nuclei having A and A — 2 nucleons. The 
trend of S2n{A) is in general rather smooth, except after having passed a major shell 
gap. There, when adding 2 more neutrons, the S2n{A) value drops significantly. This is 
an indication that a shell closure is efficient and that a spherical shell gap is existing. 
Such a sharp drop in S2niA) is observed in the Ca and Ar isotopic chains (see Figure [2]). 
On the other hand, the S and Si isotopic chains behave differently. An increase in 
S2n{A) is observed in the S chain when crossing N = 28 and in the Si chain, just before 
crossing = 28 . This indicates a clear deviation to sphericity at A^ = 28 in these 
isotopic chains. This increase of S2n{A) around A^ = 28 likely arises from the increase of 
binding energy due to correlations, which goes in concert with the onset of deformation. 
These hypotheses need to be confirmed with other nuclear properties such as the first 
2f energy and the reduced transition probability B{E2; 0^ — t- 2^). 

2.2. Evolution of E{2t) and B{E2; Of 2^^) values in the A^ = 20 - 28 nuclei 

Most generally, magic (spherical) nuclei are characterized by high E{2i) and weak 
B{E2;0i — > 2^) values. Any deviation from this trend reveals a change in nuclear 
structure. The systematics of the E{2f) and B(E2) values from '^°Ca to ^*Ca (see 
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Figure [3]) is often used to model multi-particle configurations of valence nucleons in one 
particular orbit, which is the 1^/7/2 orbit in the present case. In this framework, the 2f 
energies are maximum at the beginning of the shell (A^ = 20) where there is no valence 
neutron, and the end of the shell {N = 28) where neutron excitations can no longer 
occur in the same shell. The 2f energies are almost constant in between. In parallel 
the B(E2) values follow a bell-shape curve with minima at the two extremes of the shell 
according to the relation: 

B{E2) (X F{F -1) (1) 

where F is the fractional filling of the shell, which amounts here to F = {N — 20)/8. 

A deviation to this curve would be due to the breakdown of a spherical shell gap 
and an onset of deformation. Importantly, it could also happen when the configuration 
of the 2f state changes from a neutron origin to proton one in a given isotopic chain. 
In this case a sudden enhancement in B{E2) value can be observed as protons carry a 
much larger effective charge than neutrons (typically ep=1.5 while e„=0.5). Added to 
the E{2f) and B{E2) values, the 4^ and 0^ energies provide important complementary 
information on the evolution of the shell structure. The ratio ElAf /2f) is often used as 
an indicator of deformation, while the location of the 0^ state and the reduced electric 
monopole strength p{EO; O2 — )■ Of) are used to probe shape coexistence. 




NQutron number Neutron number 



Figure 3. Experimental E{2f) energies (left) and B{E2;0f 2f) values, also noted 
B{E2)up or B[E2)'\, (right) in the i2Mg to 2oCa isotopic chains as a function of the 
neutron number N . Data are taken from the compilation of Ref. [B], except for "^^Si [7], 
4°Si [5], "^Si 0, [10], ''^S [U, 44 Ar dU, 46 Ar [HE], and ^SAr [TS]. 

Figure [s] shows the evolution of the 2f energies and of the reduced transition 
probabilities B{E2;0f — )■ 2^) as a function of the neutron number in the i2Mg, 1481, 
igS, igAr and 2oCa isotopes. At A^ = 20, the 2+ energies are high for all isotones, except 
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for 2o^g) in which a drop in energy is observed instead. The B(E2) values are small 
and remarkably similar for all = 20 isotopes, except again for l^Mg. This picture is 
consistent with a strong shell closure N = 20 which persists from Z = 20 to Z = 14, 
and disappears below Z = 14 (see [16] for further discussions). On the other hand, 
the behavior of the = 28 isotones is different. The rise in 2f energies observed at 
Z = 16 — 18 is more modest than at = 20. In particular, the rise at ^^823 is much 
smaller than at ^^820 ■ In 148128 a drop is observed, at variance with what was observed 
in the '^^8i2o isotope. Combining this view of the 2+ energy trend and the evolution of 
S2n{A), one can surmise that the = 28 shell closure is progressively vanishing below 
Ca. The trend in B{E2) nicely confirms this statement, as discussed in the next 8ection. 

2.3. Spectroscopy of the even-Z N = 28 isotones 

In this 8ection, we discuss more precisely some properties of the even-Z A^ = 28 isotones, 
from Z = 18 to Z = 12. 

2.3.1. ^^Ar As the Ar isotopes have an open proton shell configuration, the 2^ state 
is likely to be mainly of proton origin there. The 2^ energy at the A = 20 and A = 28 
shell closures is smaller in the Ar isotopes as it comes naturally from the proton re- 
coupling inside the sd shells. However, a rise in 2+ energy is found, witnessing the 
presence of neutron shell closures. It is worth noting that two very different values of 
B{E2-0t — y 2]*^) values have been determined in ^'^Ar (Figure [3|. They arise from two 
different experimental techniques. On the one side, two experiments were carried at the 
N8CL facility to study the Coulomb excitation of ^^Ar. They led to consistent results, 
with a rather small B{E2)up value of 218(31) e^fm^ p!3]. On the other side, a fairly large 
B{E2)up value of 570 IfgQ e^fm^ was deduced from another experiment performed at the 
Legnaro facility by Mengoni et al [H]. In this work, the ^^Ar was produced by means 
of multi-nucleon transfer reaction. The lifetime of the 2^ state was determined by using 
the differential recoil distance Doppler shift method. In addition, a tentative 4^ state 
has been proposed at 3892 keV by Dombradi et al [17], leading to E(4]'"/2]^)=2.47, a 
value which lies between the limits of vibrator and rotor nuclei. 

Comparison of these two B{E2) values to models is very instructive. The 
recent relativistic Hartree-Bogoliubov model based on the DD-PCl energy density 
functional [18] reproduces very well the low value of Ref. [I3], while the 8hell- Model 
(8M) calculations of Nowacki et al [19] agree with the large value of Ref. p3]. It 
therefore seems that, though only two protons below the doubly-magic "^^Ca nucleus, 
^^Ar is still a challenging nucleus to modelize. 

Interesting to add is the fact that the B(E2) in ^^Ar |i20j is close to the lowest 
value measured in "^^Ar, i.e. at variance with the 8M predictions, while the B{E2) 
experimental result and 8M calculations agree with an enhanced collectivity for ^*Ar [20] . 
8M calculations also account for the spectroscopy of higher energy states in ^^Ar 
obtained from deep-inelastic transfer reactions, which show that its deformation is 
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likely nonaxially symmetric [15]. To conclude, a new measurement of the B(E2) in 
^^Ar is therefore called for to see if a local discrepancy between SM calculations and 
experimental results persists at = 28. 

2.3.2. ^"^S The B{E2]Qi — )■ 2^) and 2^ energy values in ^^S, determined using 
Coulomb excitation at high energy at the NSCL facility by Glasmacher et al [21], point 
to a configuration which is intermediate between a spherical and deformed nucleus. 
Evidence of shape coexistence was found at GANIL by determining the decay rates 
of the isomeric state at 1365(1) keV to the 2^ at 1329(1) keV and 0^^ ground 
state through delayed electron and 7 spectroscopy [22]. Comparisons to shell-model 
calculations point towards prolate-spherical shape coexistence. A schematic two-level 
mixing model was used to extract a weak mixing between the two configurations, the Oi" 
state having a deformed configuration (i.e. a two-neutron particle-hole 2p2h excitation 
across the A^ = 28 spherical gap) and the state, a spherical configuration (i.e. 
a zero-neutron particle-hole excitation OpO/i). New excited states in ^^S were found 
afterwards [231 121] , which lend support to the shape coexistence hypothesis. By using 
the two-proton knockout reaction from ^^Ar at intermediate beam energy at NSCL, 
four new excited states were observed [23] . The authors proposed that one of them 
has a strongly deformed configuration, due to the promotion of a single neutron {Iplh) 
across the A^ = 28 gap. It would follow that three configurations would coexist in ^^S, 
corresponding to zero-, one-, and two-neutron particle-hole excitations. The one-proton 
knockout reaction from ^^Cl gave access to other excited states and their 7 decays were 
measured using BaF2 detectors at GANIL [21]. In particular the existence of the 23" 
state at 2156(49) keV is confirmed and is likely the expected 'spherical' 2"*" state. The 
Hartree-Bogoliubov model based on the DD-PCl energy density functional [18] gives a 
relatively good agreement with experimental results, with the exception of a too large 
mixing between the 0^^ and 0^ states. The Gogny force with the DIS set of parameters 
also gives a good description of ^^S, as described in Refs. [2S1 [2S]- Shape coexistence 
has also been evidenced in the neighboring ^^S27 nucleus, which will be described in 
Section HI 

2.3.3. "^^Si As to the '^^Si28 nucleus, Bastin et al pj have established a 2+ state 
at 770(19) keV. This experiment used nucleon removal reactions from secondary beams 
centered around ^"^S at intermediate energy to produce the ^^Si nucleus and study the 7- 
rays from their de-excitation in flight. The detection of the 7-rays was achieved by arrays 
of detectors which surrounded the production target in which the reactions occurred. 
The dramatic decrease of the 2^ energy in ^^Si is a proof of the disappearance of the 
spherical N = 28 shell closure at Z = 14. This extremely low energy of 770 keV -actually 
one of the smallest among nuclei having a similar atomic mass - cannot be obtained solely 
from neutron excitations. Proton-core excitations should play an important role, which 
could in principle be evidenced by measuring the evolution of the B{E2) values in the Si 
isotopic chain while reaching A^ = 28. The B{E2) values in the i4Si isotopic chain [27] 
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Figure 4. Evolution of the B{E2) values in the Si isotopic chain. Experimental values 
are taken from Ref. [57] while theoretical values are taken from Shell-Model calculations 
using the sdpf-U interaction [19] in which monopole matrix elements contain tensor 
effects (dashed blue line) [35]. The result obtained when removing this tensor part is 
shown with a red line. The role of tensor forces is mainly seen at = 28 



seem to rise after = 20, but not as mucli as in tlie leS one (see tlie riglit part of 
Figure [3]). Wlietlier tlie B(E2) values remain small, steadily increase up to iV = 28, 
or follow a parabola cannot be judged at the present time as the quoted experimental 
error bars are too large. A reduced Z = 14 shell gap would dramatically increase the 
B{E2) values, as protons are carrying most of the effective charge in the nucleus. The 
sole decrease of the = 28 gap would barely change the 2f and B{E2) values [9]. The 
effect of the tensor force is to reduce the 2f energy and enhance the B{E2) value at 
A^ = 28 as shown in Figure Q leading to B{E2; Of 2^) ~ 430 e^fm^ in ^^Si. Without 
implementing tensor parts in the monopole terms, the B{E2) in "^^Si drops down to 150 
e^fm^. One could deduce that the studies of the 2^ state and the B{E2) value in "^^Si 
are essential to ascertain the role of tensor forces at A^ = 28. 

To conclude about the role of the tensor term in mean-field models, it requires a 
word of caution. Indeed, the B(E2) value of ^^Si has been calculated in the relativistic 
Hartree-Bogoliubov model of Li et al [18] and in the Gogny (DIS) model of Rodriguez- 
Guzman et al [25]. These two models do not contain tensor force. However their B(E2) 
values differ significantly: 200 e^fm^ in [18] and about 470 e^fm"^ in [26j. As far as the 
2^ energy is concerned, the value of Ref. [IB] is twice as large as the experimental value, 
while the one of Ref. [26] is closer to experiment. Therefore the B{E2) values should be 
measured and compared to theory in the whole Si isotopic chain to see if a significant 
increase is occurring at A^ = 28. The evolution of the sd proton orbits should be used 
as well as they influence strongly the B(E2) and E(2+) values. 

Further discussions on the implementation of tensor interaction and its role in the 
evolution of the gaps could be found for instance in Refs. [291 EQ] [3ll [32] . 
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Results of a new experimental study of the excited states of Si have just been 
published [33]. Thanks to the intense radioactive beams provided at RIKEN RIBF 
which enable 7 — 7 coincidence measurements, the most probable candidate for the 
transition from the yrast 4"*" state to the 2"*" was identified, leading to a 4^ energy 
of 2173(14) keV. Then the energy ratio, R4/2 ~ 2.9, corresponds to a well-deformed 
rotor. In addition, two other 7 lines were measured at high energy (at 2032(9) and 
2357(15) keV), which would deserve to be better characterized in order to assign other 
excited states of ''^Si. 



2.3.4- ^''Mg Some years ago, the observation of three Mg nuclei in the fragmentation 
of a primary beam of ^^Ca impinging a W target has ended the speculation about the 
location of the neutron drip-line at Z = 12 [3l]. This isotope is predicted to lie inside 
the neutron drip line in many theoretical calculations (see for instance Ref. |35] and 
references therein). The relativistic Hartree-Bogoliubov model for triaxial nuclei was 
used to calculate the binding energy map of ^°Mg in the /3 — 7 plane |T8J, which predicts 
that this very neutron-rich isotope shows a deep prolate minimum at (/3, 7) = (0.45, 0°). 
Shell model p2] as well as Gogny (DIS) [26] calculations predict an extremely prolate 
rotor as well. The identification of the first excited state of ^°Mg remains an ambitious 
challenge for the future. Conjectures about shell evolutions below ^^Si will be provided 
in Section l6l5l 



3. Evolution of sd proton orbits as a function of the neutron number 

3.1. The 7t2s 1/2 and it Id^/ 2 orbits 

The change of structural behavior between the iV = 20 and = 28 isotones can be 
partly ascribed to the evolution of the proton Single-Particle Energies SPE [3^. Using 
the pick-up reaction {d,^B.e) from stable 2oCa targets, the evolution of the 7r2si/2 and 
7r2(i3/2 spacing has been revealed in the 19K isotopic chain. The filling of eight neutrons 
in the 1^2/7/2 orbital induces an enhanced binding energy of the 77^3/2 orbit as compared 
to the 7r2si/2 one. The spacing 7r2si/2 — 7rlc?3/2 derived from SPE values drops from 
2.52 MeV to -0.29 MeV, i.e. the orbits have crossed at = 28. This is likely to be 
due to the fact that the monopole interaction 1^1^3/21/7/9 1 more attractive than the 



2 1 



l^2si/2i/7/2 1 ^'^^^ derived that 

J-a3/2l/7/2 ^^1/2-1-/7/2 \ ' 

The fact that ^1^3/31/7/2 significantly more attractive than ^251/21/7/2 could 
quahtatively be ascribed to the fact that the proton 1^3/2 (^=2, £ I) and neutron I/7/2 
(£=3, i t) wave functions have the same number of nodes and have attractive tensor 
terms as the proton and neutron spins are anti-aligned. On the other hand for ^251/21/7/2' 
the numbers of nodes differ and a large difference in the orbital momentum value is 
present between the two wave functions , making this monopole term weaker. Taking 
into account the monopole matrix elements solely, the evolution of [E{l/2~^) — E{3/2~^)] 
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Figure 5. Calculated (black dashed line) and experimental (red line) values of 
[£'(1/2+) - £(3/2+)] along the K, CI and P chains (adapted from Ref. |Mj)- The 
fixed monopole-driven trend (blue line) given by Equation [S] accounts well for the 
global evolution of [£(1/2+) — £(3/2]^)], while correlations are required at mid-shell. 
Note that for the P chain the ordering of the 1/2+ and 3/2+ states has changed as 
they correspond to particle states while they are hole states in the CI and K chains. 



between = 20 and = 28 would be linear as a function of the number of neutrons 
X leading to: 

[E(l/2+)-E(3/2+)]2o+. = [i?(l/2+)-E(3/2+)]2o+x(\/2;^^^/^-l^^;^^^^^J(3) 

The experimental evolution of [E{l/2'i) — E{3/2i)] deviates at mid shell from the linear 
monopole trend shown in Figure |5| This deviation is due to pairing and quadrupole 
correlations which already engage at = 22, as soon as the single-particle states vrsi/2 
and 7id3/2 come close enough to each other. Shell model calculations using the sdpf 
interaction [37] well reproduce these correlations [381 [39], ^ shown by the black dashed 
line in Figure |5} 

The strong reduction between the proton 3/2+ and 1/2+ states is also found around 
= 28 in the 17CI and 15P isotopic chains. This demonstrates that the change 
in monopole interaction plays a decisive role in bringing the 7rsi/2 and 71^3/2 states 
degenerate at A^ = 28. This has a profound consequence on the evolution of collectivity 
between A^ = 20, where a sub-shell gap exists between ttsi/2 and 71^3/2, and A^ = 28, 
where this sub-shell has vanished completely. As discussed above for the K isotopes, the 
evolution of the experimental 3/2^ and 1/2^^ states for the CI and P isotopes is distorted 
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Figure 6. Left: Evolution of the proton dr,/2, Si/2 and orbits in the K isotopic 
chain derived from experimental data. Right: Results of HFB calculations using the 
Gogny DIS interaction (see text). 

by pairing and quadrupole correlations (see Figure [s]), which are also well accounted for 
by shell- model calculations 

3.2. The 71(15/2 orbit 

As shown before, the evolution of the Z = 14 shell gap is crucial for providing enhanced 
correlations in the N = 28 nuclei far from stability, as well as to probe the effect of 
tensor forces. However, as the single-particle strength is significantly spread in the ^^P 
nucleus, the size of the Z = 14 gap can hardly be extracted there from experimental 
data. In order to extract more accurately the change in the size of the Z = 14 gap 
when filling the neutron fjf2 orbit from = 20 to = 28 we could in principle look 
at the evolution of the si/2 and 6/5/2 single-particle energies in the K isotopic chain as 
deformation is not developing there. 

However, one additional problem arises from the fact that in the K chain the ^5/2 
orbit appears to be more bound by about 2.5 MeV than in the P isotopic chain. It 
follows that it is very hard to study the evolution of the ^5/2 single particle energy in 
the K chain as the ^5/2 strength is spread into many states which carry a small fraction 
of it (see the discussion in Section 4.2.3 of Ref. |3] and its Tables 2 and 3). With this 
important word of caution in mind, we can nevertheless discuss the evolution of the 
binding energies of the 7r(i5/2, vrsi/2 and vr(i3/2 orbits between ^^K2o and ^''K28 in a semi 
quantitative way. The experimental part of Figure |6] displays a reduction of the proton 
d5/2 — splitting by about 1.7 MeV between A^ = 20 and A^ = 28. When adding 
neutrons in the /7/2 shell, the two d orbits become more bound. This is mainly due to 
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the attractive proton-neutron interactions. The fact that the gain in binding energy is 
larger for the 1^3/2 orbit than for the 1^5/2 one comes from the fact that the ^1^5/21/7/2 
monopole is weaker than ^1^3/21/7/2- likely comes (at least partly) from the fact that 
^^3/2 1/7/2 ''^Wr/2 1/7/2) J^oiiopols contains an attractive (repulsive) tensor part as the 
proton and the neutron have anti-aligned (aligned) spin orientations. 

Note that the size of the Z = 14 gap, which is almost similar for the two isotones 
■^^K and ^^P (the properties of ^^P are discussed in Section 3.4), is an essential ingredient 
for providing a similar structural behavior at = 20 from ^°Ca to ^^Si. 



3. 3. Theoretical predictions of the evolution of the three orbits in the K isotopic chain 

Interesting is to compare this variation of proton orbits to calculations obtained with 
Hartree-Fock (HF) self-consistent calculations using the DIS Gogny force [lOll^ . The 
binding energies of the single-particle levels of ^"Ca and '^^Ca have been obtained 
from HF calculations. As usual, the HF equations are firstly solved by iterative 
diagonalization of the HF Hamiltonian. Then the single-particle energies are defined as 
the eigenvalues of the self-consistent one-body Hamiltonian, obtained after convergence. 
In this approach, a change of the Si/2 — splitting is found from = 20 to = 28, 
but twice as small as experimentally. No change of the ^5/2 — ^3/2 sphtting is obtained 
(see the right part of Figure |6]) as tensor interaction is not included. Moreover, as 
these orbits are well bound no reduction of the spin-orbit interaction associated with 
an extended surface density profile is foreseen. We conclude that a further reduction of 
the d^i2 — (^3/2 splitting is required to match the experimental results. 

3.4. Comparison of^^P2o and'^^P28 

The development of collectivity in the igS and uSi isotopic chains depends on the amount 
of proton excitations across the Z = 14 gap, from the deeply bound proton ^5/2 to the 
degenerate ^3/2 and si/2 orbits (see the bottom part of Figure [T]). When possible, the 
■5i/2-'^5/2 cross-shell excitations naturally bring quadrupole excitations. The comparative 
studies of the ^^P2o and ^^P28 nuclei gives valuable information on the reduction of the 
Z = 14 gap and the increasing role of correlations across it. 

The energy of the three orbits, 7id5/2, vrsi/2 and 'n'd3/2, and their occupancy values, 
(2 J + 1) X SF (where SF is the spectroscopic factor), have been obtained in ^^P2o by 
Khan et al |12] by means of the '^^S2o('^,^He)^^P2o reaction. The major results are 
shown in the top part of Figure [7| The energy spectrum of ^^P2o exhibits few levels, 
with rather large spacing. The sum of the occupancies of the two first states Si/2 and 
(^3/2 is about 2, while most of the occupancy is in the Si/2 orbit. As the Si/2 — ^3/2 
spacing is not large enough (about 2.4 MeV) to smother pair scattering, some leak of 
occupancy from the Si/2 to the upper ^3/2 orbit is present. The full 0^5/2 strength is 
shared in three major fragments between 3.8 and 5.2 MeV excitation energy, leading to 
a mean value of the Z = 14 gap of 4.6 MeV. The picture is there rather simple, with 
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Figure 7. Top: Occupancy values of the proton ^5/2^ Si/2 and ^3/2 orbits are given for 
'^^P2o- They were derived from spectroscopic factors obtained in the ■^®S2o(d,'^He)^^P2o 
reaction [32] • The first excited states of ^^P2o are shown in the right part of the figure. 
Bottom (adapted from Ref. [33]): calculated occupancy values for '*'^P28 derived from 
the '^''828 (-Ip) ^'^P28 reaction. The level scheme of '*'^P28 shown in the right part of 
the figure contains many 5/2"*" states, witnessing a large spreading of the proton d^ji 
strength in this nucleus. 



a sequence of ^5/2, Si/2 and ^3/2 orbits bound by about 16.8 MeV, 12.2 and 9.8 MeV, 
respectively. 

The situation gets more complex in the ^^P28 nucleus, which was studied by means 
of one proton knock-out reaction ^''828 (-Ip) "^^1^28 by Riley et al [33]. At first glance 
it is clear in Figure g that the energy spectrum of P28 is more compressed and 
contains many more levels than ^^P2o does. This compression of levels is relatively 
well reproduced by shell-model calculations using the interaction of Utsuno et al [H] . 
If the sum of the si/2 and c/3/2 occupancies still leads to about 2, the sharing of occupancy 
reflects the quasi degeneracy of these two orbits: since the (^3/2 orbits contains twice as 
many sub-states as the si/2 does its occupancy is implicitly twice as large. As to the ^5/2 
strength, it is spread over many states among which only the ones below 2 MeV have 
been observed experimentally. Noteworthy is the fact that the ^5/2 strength already 
starts at low energy and a significant fraction is already calculated below 3 MeV. This 
high density of states at low energy as well as the spreading of the occupancy value 
are at variance with the picture observed in ^^P. In order to obtain a significant ^5/2 
strength below 3.5 MeV in ^^P, the spherical Z = lA gap has to be slightly reduced 
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between = 20 and = 28 [H E] . 
3.5. Summary 

The development of collectivity in the neutron-rich = 28 isotones is partly due to 
the reduction of the spacings between the protons ^5/2, Si/2 and 6/3/2 orbits as soon as 
the neutron fY/2 orbit is completely filled. As these orbits are separated by 2 units of 
angular momentum, quadrupole (E2) collectivity is naturally favored. Added to this, a 
reduction of the N = 28 shell gap would reinforce this tendency to deform. Mean-field 
and shell-model theories agree on this description. We note that from the evolution 
of the proton sd orbits, mean-field theories need an additional spin-dependent term to 
reduce the ^5/2 — (^3/2 splitting in order to match the experimental value. 

It is important to note that a systematic study |l5l HH] of the low-lying structures 
of many isotopic and isotonic chains was performed within the shell model approach, 
using the interaction of Nowacki and Poves [I9]. In these works, a particular focus is 
also made on the cases with unpaired nucleons. Indeed the properties of the N = 27 
and A^ = 29 isotones give a deeper understanding of the evolution of the A^ = 28 shell 
gap, as detailed in the two next Sections. 

4. Structural evolution viewed from the N = 27 isotones. 

The possible development of collectivity can be looked at in the N = 27 isotones, by 
comparing the characteristics of their first states. Shell-model calculations [lEj have 
been performed by Gaudefroy et al using the interaction of Nowacki and Poves [19]. 
The results are drawn in Figure [8] In this plot, the configurations of the ground and 
the first excited states are displayed using a bi-dimentional representation, where the 
proton (neutron) configuration is shown on the y (x) axis. A^^^ {N^) gives the number of 
proton (neutron) excitations above the OpOh (Oplh) core configuration and the size of 
the squares in the N-j^ — Ny representation gives the intensity of each component. 

About 90% of the g.s. configuration of the spherical ^^Ca nucleus corresponds to a 
neutron hole inside the /7/2 shell (Opl/i, Ny=Q) and a closed proton core (OpO/i, Nt,=Q). 
The first-excited state with = 3/2" at 2.02 MeV is expected to involve the promotion 
of one neutron in the upper P3/2 shell with two neutron holes coupled in the /7/2 shell 
{lp2h, Ni,=l). This is in agreement with the theoretical results showing that the first- 
excited state at 2.07 MeV corresponds mainly to the pure excitation of one neutron 
{Ny=l) and to a minor extent, of three neutrons (A^y=3). This state, with a closed 
proton configuration, would be found well above the ground state of all the N = 27 
isotones, except if the A^ = 28 shell gap is reduced and correlations dominate. 

The configurations of the g.s. and excited states in ^^Ar have been studied using 
^^Ar(d,p)^^Ar |15] and '^^Ar(-ln)^^Ar [17] transfer and knock out reactions, respectively. 
Very good agreement is found between the experimental results, excitation energies and 
spectroscopic factors, and the theoretical predictions [IS] . The drawings of Figure M 
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Figure 8. Squared wave functions of tlie first 7/2~ and 3/2~ states in the N ~ 27 
isotones represented in the proton (Nt^) versus neutron (iV^) particle-hole configuration 
plane (see text for details), adapted from Ref. ^5]. Some of the excitation energies are 
given in the text. 



show that as compared to ^^Ca into which proton excitations are hampered by the 
Z = 20 gap, proton excitations are naturally present in ^^Ar inside the Si/2'^3/2 proton 
states which are degenerate in energy at ~ 28, as discussed in Section |3.1[ These 
proton excitations occur as Iplh {N^^ = 1) or 2p2h (A^^ = 2). They both contribute 
to the largest fraction of the B[E2) values in the Z = 16, 18 nuclei. In line with these 
features the ^^Ar ground state is still dominated by a configuration similar to the one 
of ^^Ca, but correlations leads to a more mixed wave function. The state corresponding 
to the first-excited state in ^''Ca is more likely the second-excited state in ^^Ar (3/2^) 
observed at 1.42 MeV and predicted at 1.22 MeV, for which the N^=0 component does 
not exist either. The configuration of the (3/2|f) state at 550 keV is in between the 
7/2~^ and the 3/22 configuration, which already witnesses the sharing of the p strength 
among several states. 

Below ^^Ar, an inversion between the 7/2~ and 3/2~ states is predicted (see 
Figure [s]). An isomer has been discovered in ^^S by Sarazin et al at 320 keV ^5] . 
It was interpreted with the help of SM calculation to have a 7/2~ spin value and to 
decay to the 3/2~ ground state by a delayed transition. More recently, the g factor 
of this isomer was measured [19] establishing its spin-parity values, I'^ = 7/2^, and 
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implying a rather spherical shape for this state. Nevertheless, the value of its quadrupole 
moment is larger than expected for a single-particle state [50], implying that it contains 
non-negligible correlations that drive the state away from a purely spherical shape, in 
agreement with SM calculations (see Figure [s]). In addition, an intermediate-energy 
single-neutron knockout reaction was used to characterize other excited states of ^^S 
|51j . Two of them are proposed to be members of the rotational band built on the 
deformed 3/2^ ground state, strengthening the case for shape coexistence in ^^S. 

Note that this inversion of the 7/2~ and 3/2^ states, and the spreading of the 
wave function among neutron and proton excitations is, according to SM calculations, 
persisting in the ^^Si27 isotope. The spectroscopy of ^^'^^Si has been carried out by Sohler 
et al |52] using the in-beam 7-ray spectroscopy method from few nucleon knockout 
reactions. The observation of low lying states in ^^Si call for a reduction of the N = 28 
and Z = 14 shell gaps which induce the lowering of the intruder neutron 3/2~ state 
while going from Z = 20 to Z = 14. The energy of the only 7 line, 672(14) keV, 
observed in ^^Si is significantly lower than the one of the first excited state in ^''Ca 
(2014 keV) suggesting a deformed ground state for ^^Si. However, as compared to shell- 
model predictions, some states expected at low energy have not been observed in ^^Si. It 
was suggested in Ref. [52] that a low-energy isomer would be present there. As isomers 
could not be evidenced with in-beam spectroscopy, another technique has to be used to 
reveal its existence. The search for this isomer and/or low-energy states would bring 
important pieces of information to confront to SM description. In Figure |8] one can 
observe that for "^^Si many squares have equal size, pointing to large mixing of proton 
and neutron configurations. 

As mentioned above, the inversion between the natural (2^/7/2)^^ and the intruder 
{^P3/2)~^^ configurations of the ground states of the N = 27 isotones occurs between ^^Ar 
and ^'^S. The properties of the ground state of ^^Cl therefore allow us to delineate the 
exact location of the inversion. This odd-odd nucleus has been populated by using the 
single neutron-knockout reaction from ^^Cl28 at intermediate beam energies at MSU [53] . 
The momentum distribution of the recoiling ^*^C1 nuclei was analyzed for the direct 
population of its ground state, the best fit being obtained for the removal of an £ = 1 
neutron. This means that (i) there is a significant (1^^3/2)^^ component in the ground 
state of ^^Cl and (ii) this neutron orbit is involved in the ground state of ^^Cl. Therefore 
the intrusion of the i'P3/2 orbit from above the N = 28 shell closure occurs already for 
Z = 17. This is confirmed by the value of the g factor of the ground state of ^^Cl. 
This measurement, done at the LISE fragment separator at Ganil using the /3 nuclear 
magnetic resonance technique |54j, indicates that (^(^^Cl) is significantly lower than 
g{'^^K). The two g values are well reproduced by SM calculations. While the main 
configuration of the ground state of ^^K is dominated by the spherical (vr(i3/2)^(i^/7/2)^ 
configuration, the wave function of the ground state of ''^Cl is very fragmented, its 
second most-intense component, (7rsi/2)^(^'P3/2)^, explaining why the g factor of "^^Cl is 
reduced 
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5. Structural evolution viewed from the = 29 isotones. 

In the previous sections, the erosion of the = 28 shell gap below the doubly magic 
2QCa has been probed through properties of atomic masses, energies of excited states 
and reduced transition probabilities, as well as through the spectroscopy of = 27 
isotones. The most direct way to evidence this erosion is to determine the evolution of 
the single particle energies between the 20^^ is^^ isotones. The study of '^'^At29 
has been carried out by the ^gAr((i,p) transfer reaction in inverse kinematics at the 
GANIL/SPIRALl facility [55j. The obtained results were used to investigate the change 
of the A^ = 28 gap from 2oCa to isAr. From the Q value of the transfer reaction, the 
A^ = 28 gap was found to be of 4.47(9) MeV in "^^Ar, which is 330(90) keV smaller than 
in ''^Ca. Transfer to excited states were used to determine the energies and spectroscopic 
factors of the neutron ^3/2, pi/2 and /5/2 states in f|Ar. The fact that only part (about 
3/4) of the strength for these states has been observed indicates that some correlations 
takes place already in ^^Ar, after the removal of only 2 protons from the doubly magic 
^*Ca nucleus. In particular it was found in Ref. [55] that the P3/2 orbit is already 
partially occupied through Ip — Ih excitation across the A^ = 28 gap, from the /7/2 
orbit. The agreement between experiment results and SM calculations was rather good 
using the monopole terms of the sdpf |37] interaction. They were however slightly 
adjusted to match existing experimental data in this mass region [56]. By applying this 
technique, the full particle strength of the neutron /7/2, P3/2 and pi/2 and /5/2 orbits 
has been determined in ^^Ar. The resulting SPE are compared to those of the 20^^ 
isotone [STJ |58] in Figure [9] It is found that from ^''Ar to ^^Ca the orbits in which the 
angular momentum is aligned with the intrinsic spin such as /7/2 and P3/2, become 
relatively more bound than the /5/2 and pi/2 orbits where the angular momentum and 
intrinsic spin are anti-aligned (ii). In particular, one can see in Figure |9] the asymmetry 
in gain of binding energy between the /7/2 and /5/2 orbits. 

Bearing in mind that the ds/2 and Si/2 orbitals are quasi-degenerate at A^ = 28, (see 



Section 3.1 ), the addition of 2 protons between ^^Ar and ^^Ca occurs in an equiprobable 
mannei[§]in these orbits (i.e 1.33 in 1^3/2 and 0.66 in 2si/2). Therefore modifications of 
neutron SPE arise from proton-neutron interactions involving these two orbits [HS] and 
the change of the A^ = 28 shell gap {SG) can be approximated to : 

5G-(28) ^ 1.33(^-^2.3/. - ^1^.1/./.) + 0-66(^2^.;.2p3/. - ^i^;,!/,/,) (4) 
Similarly, changes in the p and / SO splitting express as: 

6S0 = 1.33(1/,- - VrZ^^,^) + 0.66(1/- - ,,,) (5) 

As regards the / states, the experimental change in SO splitting 6S0{f) was ascribed 
in [55] to the fact that the monopoles contain attractive and repulsive spin-dependent 
terms (which could be tensor terms) : 

6S0{f) ^ 1.33(1/,- - V^-;^,,^^J (6) 

§ When taking into account the correlations expected in non doubly-magic nuclei, the occupation rates 
of the orbits are slightly modified, see for instance the numbers given in Ref. [5 6) 
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Figure 9. Left: Neutron single-particle energies (SPE) of the fp orbitals for the 
|gAr29 and 2oCa29 nuclei. Right: Schematic view of the proton-neutron spin-dependent 
interactions involved to change the / (top) and p (bottom) SO splittings derived from 
the experimental data on "^^Ar in Refs. [551 156) . 
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As shown in Figure 9 the spin-dependent parts of the monopoles V amount to V{^, 



3/2 1/7/2 

=-210 keV and Vf^^^^ ]^j^^^=-|-280 keV, respectively. They amount to about 20% of the 
total monopole term The change of the p SO splitting 6S0{p) was principally 

assigned in [SB] to the removal of a certain fraction of 2si/2 protons which depletes the 
central densit}|ji]of the nucleus: 



6S0{p) ~ 0M{V11 



/22pi/2 



ypn ^ 



(7) 



As shown in Figure |9j the corresponding spin-dependent part of the monopole terms 
was extracted to be of -1-170 keV and -85 keV for V^^"^ ^^^^ and ^^"^^ 2^3/2' respectively. 

The variation of the single-particle energies in terms of monopole interactions could 
be pursued towards the ^4Si nucleus in which, as compared to 2oCa, about 4 protons 
have been removed in the ^3/2 and 2 in the Si/2 orbits. We note here that protons are 
partly removed in the orbit as well. Using the monopole matrix elements derived in 
the interaction, changes in the spacing of the fp orbits is foreseen when the ^3/2 and si/2 
orbits are empty. This leads to reductions of (i) the = 28 gap by 0.33 x 3 ~ 1 MeV, 
(ii) the / SO splitting by ~ 4(0.28 + 0.21) ~ 2 MeV and (iii) the p SO splitting 
by ^ 2(0.17 + 0.085) ~ 0.5 MeV. These global reductions are expected to reinforce 

II Having an angular momentum £=0, the Si/2 orbit is peaked in the center of the nucleus. When no 
proton occupies the Si/2 orbit, the central density is therefore depleted as compared to nuclei in which 
this orbit is filled. This property will be used in Section [63| to study the density dependence of the SO 
interaction. 
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particle-hole excitations across = 28, which are of quadrupole nature. 

It is important to note that, in the shell- model approach, the monopole terms do 
not vary when approaching the drip-line. Therefore when reaching the continuum the 
present approach and the concepts that have been used will be inappropriate. The 
proximity of continuum could already modify the binding energy of the /5/2 orbit in 
^^Ar, hereby possibly altering the role of tensor interaction derived from the present 
experimental data. The role of continuum should be quantified in this region. 

6. Questions and perspectives 

6.1. Introduction 

As described in the previous sections, a wealth of spectroscopic information has been 
obtained from the study of the A?" ~ 28 nuclei from Z = 20 down to Z = 12. Several 
important questions related to nuclear forces arose from these studies. In the following 
we propose to address some of them: 

• What is the role of the three-body forces in providing the SO shell gaps ? Is there 
a general picture that is present throughout the chart of nuclides ? 

• Can we disentangle the respective roles of the tensor and the spin-orbit components 
to account for the disappearance of A'^ = 28 far from stability? Can we better 
understand the spin-orbit force ? 

• Would other SO magic numbers (14, 50, 82 and 126) be disappearing far from 
stability as = 28 does and if not why ? 

• Which forces come into play in nuclei located below ^^Si, ^^Ni and ^^^Sn? 

• Are proton-neutron forces significantly modified when approaching the drip-line ? 

We propose to discuss these questions using rather qualitative and phenomenological 
arguments based on experiments results. 

6.2. What is the role of the three-body forces in providing the SO shell gaps ? 

6.2.1. Introduction In this Section we investigate the role of neutron-neutron 
interactions to create the neutron SO shell gaps. The three-body forces are needed 
to create these gaps, which would not be obtained if realistic two-body forces would be 
used exclusively. We start with the study of the A^ = 28 shell gap and generalize this 
study to other SO shell gaps as A^ = 14, A^ = 50 and A^ = 90. We conclude that the 
three-body forces seem to play a similar and essential role to create the major shell gaps 
that leads to SO magic nuclei. We propose an empirical rule to predict the evolution of 
high-j orbits in superheavy nuclei as well as to create a sub-shell gap N — 90, which 
could be of importance for the r-process nucleosynthesis. 
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Figure 10. Evolution of neutron binding energies in the O (left), Ca (middle) and Ni 
(right) isotopic series. In each chain, the value of the neutron shell gap grows when 
adding neutrons in the high-j orbit (^5/2, /7/2: 99/2 respectively). In the cases of the O 
and Ca chains, experimental values are taken at the two extremes of the binding energy 
values and a straight line is drawn in between, according to the monopole-driven trend. 
As for the Ni chain, details are found in the text. 



6.2.2. Neutron- neutron forces to create the = 28 shell gap in Ca Search for 
fp single-particle have been carried out using transfer reactions [59^ i58j at = 20 and 
iV = 28 by studying the particle and hole strengths around the ^°Ca and ^^Ca nuclei, 
respectively. These experiments show that the = 28 shell gap strongly increases 
when adding neutrons into the fj/2 orbit, as the neutron fj/2 (P3/2) orbit gains (looses) 
binding energy. The gain in energy of the = 28 gap from A^ = 20 to A^ = 28 is about 



(5G'"'"(28) = 3 MeV (see the middle of Figure 10 ). The mechanism which creates the large 



N = 28 shell gap, between the /7/2 and P3/2 orbits, is likely due to a strongly attractive 
^7/2 /7/2 ^ repulsive Vj^^^ ^^^^ monopole termj^ The increase of the A^ = 28 shell 
gap, (5G'™(28), due to the neutron-neutron interactions can be written as a function of 
X (the number of neutrons added in the /7/2 orbit) as : 

^G-m=m:i.vj-^^-m;.f.,.\ (8) 

We note that the (x — 1) term applies when neutrons occupy the same shell. 

Effective two-body interactions derived from realistic interactions could not account 
for this increase of the A^ = 28 gap. For instance, it has been recently shown by Holt et 
al [UD] that the A^ = 28 gap is almost constant using two-body forces and grows when 
three-body forces are taken into account. A more recent work which uses interactions 
from chiral effective field theory, leads to the same conclusion: the three-nucleon forces 

\ We note here that the Y™ are effective monopole terms. Even though they are written in a two-body 
form, they may contain three-body components. 
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create a = 28 magic number [6T] and generate the shell closure in ^*Ca. More 
generally it was proposed by Zuker |62] that the implementation of a three-body force 
can solve several deficiencies in nuclear structure. 

6.2.3. Neutron-neutron forces to create the = 14 sub-shell gap in O The 

= 14 subshell gap, which is located between the ^5/2 and Si/2 orbits, also comes 
from the SO splitting and shares similar properties as the = 28 one. As shown in 
Figure 10, it grows by about (5(7"" (14) =2. 7 MeV from -'^'^O to ^^O as the neutron ^5/2 



orbit is filled. Similar to the A^ = 28 gap. Equation [8] applies and the large value of the 
A^ = 14 gap at A^ 
monopoles terms. 



A^ = 14 gap at A^ = 14 is due to the strongly attractive (repulsive) VJ^J^^ ^^^^ si/2 



6.2.4. Could we predict the size of the A^ = 50 shell gap in Ni? From a 
phenomenological point of view, we can reasonably expect that a similar increase in 
binding energy of major shells having £ and s values aligned i c?5/25 f7/2, 99/2, 

hu/2, will be occurring during their filling. This generally leads to an increase of shell 
gap by: 

sG--{so) = x[v,;^,_,)^] - (x - i)[vi;i] (9) 

In particular, we could look at what is expected for the A^ = 50 shell gap from ^^Ni 
to ^^Ni, where fewer experimental data are available so far. Based on Equation [o] the 
gap formed between the gg/2 and ^5/2 orbits is expected to grow as the l5'9/2 orbit is 
filled. The relevant effective monopole terms {V^JJ^ g^^^ ~ -200 keV and V^""^ ^^^^ ^+130 
keV) can be extracted from the spectroscopy of s^'^^Zr (see Figure 6 of Ref. |[3j). When 
renormalizing the monopole terms to A ^ 75 (using the hypothesis that V scales with 
A^^^^) we find that the A^ = 50 gap should increase by about 3.1 MeV between A^ = 40 



and A^ = 50. The right-hand side of Figure [10] displays the tentative change of the 
A^ = 50 gap in the Ni isotopic chain. The Sn value of ^^Ni as well as the preliminary 
energy of the 2^5/2 state at ~ 2.6 MeV derived from Ref. [HS] are used. Note that a 
correction in binding energy is applied to ^n(^^Ni). Indeed it is likely that the ^^m^o 
ground state has a mixed contribution with the 02" state at 1.77 MeV. Taking a mixing 
of 50%(i/5'9/2)^ + 50%(i/pi/2)^ it implies that before mixing the (pi/2)^ configuration 
of ^^Ni is likely less bound by about 0.9 MeV. Thus starting from 2.6 MeV in ^^Ni, the 
A^ = 50 shell gap would reach (2.6+3.1=) 5.8 MeV in '^^Ni. This is somehow larger 
than what extrapolated from the A^ = 50 trend of binding energies [M] and with the 
SM predictions of Ref. ^3]. Note that a small change in one of the V"'^ values given 
above can drastically impact the value of the increased gap, as 10 neutrons [x =10) are 
involved here between A^ = 40 and A^ = 50 (see Equation [9]). 

If this later increase of shell gap is confirmed, this will remarkably show that similar 
forces (including the three-body term) are at play in the various regions of the chart 

This hypothesis comes from the comparison with experimental results obtained for 4QZr, showing that 
the configurations of its ground state and its 0^ state at 1.76 MeV are 50%(7r(7g/2)^ + 50%(7rpi/2)^- 



The N = 28 shell closure: a test bench for nuclear forces 



22 



of nuclides to produce such large SO shell gaps. Even more a priori surprising is the 
increase of all the = 14, = 28 and = 50 shell gaps by a similar value of about 
3 MeV. This comes from the fact that two competing terms are involved to enlarge the 
gap in Equation [9} The two monopoles decrease with A~^^^ as the nucleus grows in size, 
but they are multiplied by the number of neutrons x in the i orbit, which increases 
between = 14 (x = 6 in 4/2), A^ = 28 (a; = 8 in /y/s), and A^ = 50 (x = 10 in 
99/2)- We could tentatively propose as an empirical rule that, during the filling of the 
i t orbit, the variation of the major shell gaps of SO origin amounts to about 3 MeV. 
Though rather crude, this empirical rule could be used to predict the evolution of the 
(j t? (j " 2) f) orbits as a function of the filling of the j ^ orbit. However, the pairing 
effect will dilute the orbital occupancy in heavier systems, as mentioned below for the 
A = 132 region of mass. 

6.2.5. Evolution of the A^ = 82 shell gap in Sn As for the next shell A^ = 82, 
we probably expect the same increase of the /iii/2 — /7/2 spacing while filling the /iii/2 
orbit. However this increase of gap could hardly be evidenced experimentally as the /iii/2 
orbit is located close to two low-£ orbits, si/2 and ^3/2. So the pairing correlations dilute 
the occupancy of the high-j orbit with the two others, implying that the /iii/2 orbit 
is gradually filling when adding 18 neutrons from A^ = 65 to A^ = 82. Consequently 
the corresponding evolution of the /iii/2 binding energy is not only governed by the 
attractive Vi"" i, interaction, but also by the repulsive ones, V^"- „ and Vi"" ^ . 

'^11/2 '^11/2 ' J f 1 ftii/2 Sl/2 "11/2 "3/2 

Therefore, in this case, the SO empirical rule could be used as a guidance to constraint 
the spacing between the i t ^md {i — 2) orbits at the mean-field level, before taking 
the role of correlations into account. It is important to point out that none of the 
self-consistent calculations reproduce the experimental location of the i^/iii/2 orbit in 
the doubly magic ^^^Sn, i.e., when this orbit is completely filled. All mean-field models 
using Skyrme and Gogny forces or relativistic NL3 and NL-Z2 forces predict that this 
high-j orbit is located just below the N = 82 gap, the 1^^3/2 orbit being more bound 
by about 1 MeV (see Figure 7 of Ref. [66]), while the ground state of ^^^Sn has spin 
3/2"*" and its first excited state at 65 keV has spin ll/2~. One surmises that the lack of 
three-body forces is responsible for the discrepancy. Indeed these forces would strongly 
enhance the binding energy of the filled Z//111/2 orbit. 

6.2.6. A new sub-shell closure at A^ = 90 ? An interesting effect of the three-body 
forces would be the creation of a gap at ^''^Sngo, between the 2/7/2 and 2p3/2 shells, 
as proposed by Sarkar et al [67j. Indeed the neutron-neutron monopole terms (and 
possibly the three-body terms) which intervene here are somehow similar to the ones 
at A^ = 28. The only modifications are that the nodes of the wave functions differ, and 
that the monopole term for A ~ 48 has to be downscaled to account for the reduction of 
the interaction energy as the nucleus grows in size at A ~ 140. When using Equation [8} 
the suitable downscaled monopole terms, and the known 2/7/2 — 2^3/2 spacing of 0.85 
MeV in ^^^Sn [HS], ES] it is found that the A^ = 90 gap is expected to grow by about 2.2 
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MeV while filling the 2/7/2 orbit. This leads to a subshell gap of about 3.1 MeV at ^ Sn. 
If established, this subshell gap would bring an additional credit to the mechanism of 
shell gap creation by three-body forces in another region of the chart of nuclides. In 
addition, this subshell could bring a local waiting point around ^^°Sn in the r-process 
nucleosynthesis. 

In summary, the three-body forces are needed to create the neutron SO gaps. While 
these gaps are not obtained when using microscopic effective interaction (see for instance 
the discussion recently done by Smirnova et al [70]), the phenomenologically adjusted 
values of the TBME or the explicit implementation of the three-body forces in SM 
approaches generate the expected neutron SO gaps. 



6.3. Could we disentangle tensor and spin-orbit forces in the N = 28 region ? 

6.3.1. Introduction The study of the nuclear structure around = 28 has revealed, 
using the shell model approach, that the tensor and two-body spin orbit interactions 
are both at play to reduce the proton and neutron gaps far from stability. As proton 
orbits are degenerate at = 28, these two components could not be disentangled. To 
separate these components, we propose in this Section to study the evolution of p SO 
splitting between ^gS and '^4Si (a bubble nucleus). We find a sizable reduction of the 
SO splitting there, which is ascribed to the spin-orbit interaction (and not to tensor 
or central forces). This assumption comes from the fact that between these two nuclei 
protons are mainly removed from the 2si/2 {d- = 0) orbit, which does not exhibit a 
specific orientation between the orbital momentum and the intrinsic spin value. 

In a second part we propose to use the bubble nucleus ffSi to test the density 
and the isospin-dependent parts of the spin-orbit interaction in mean-field approaches 
which have never been tested so far. These properties of the spin-orbit interaction are 
of minor importance in the valley of stability but are crucial at the neutron drip-line 
and in superheavy nuclei. 



6.3.2. Studying the two-body spin-orbit interaction in the SM approach As 

mentioned in Section 3.1, the proton si/2 and 6/3/2 orbits are degenerate at A^ = 28. 
The reduction of the p SO splitting at A^ = 28, between ^^Ca and ^''Ar, is due to the 
combined effects of proton-neutron interactions involving several monopole terms having 
potentially tensor and spin-orbit components. Applying Equation [5] to the p orbits, the 
SO splitting is: 

SSOip) = 1.33(F- ^2.,,, - ^r3/.2p3/J + 0-66(^2^.;,2p,/, - Vll^,,^^J (10) 

The two contributions of this Equation can be estimated from the study of the evolution 
of the p splitting in A^ = 21 nuclei. As the Si/2 and ^3/2 proton orbits are separated by 
about 2.5 MeV (see ^^P in Figure [s]), one can therefore assume that these proton orbits 
are filled sequentially in the A^ ~ 20 isotones. The pairing is providing a small dilution 
of occupancies among the two orbits as shown in Ref. |12] and Figure [Tj Between 2oCa2i 
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and feS2i, four protons are removed. The evolution of the [z/2p3/2 — z/2pi/2] SO sphtting 
writes: 

5SO{p)..4{V,X,,,,,,^-Vri^,^J (11) 
By fooking first at the evolution of the neutron [2^3/2—2^1/2] SO splitting between '^^Ca2i 



and ^''S2i in Figure 11 , one investigates the role of proton-neutron interactions involving 
the 1(^3/2 protons. It is seen that the p SO splitting amounts to 2 MeV in the two nuclei. 
It follows that the SO splitting does not change between "^^Ca and ^^S while 4 protons 
are removed from the shell, leading to SSO{p) — and ^id^^2^Pi/2 ~ ^^d,^/2'^Ps/2' 

We 

note that taking the full observed p strength, the conclusion does not change. Indeed 
the SO splitting amounts to about 1.7 MeV in the two nuclei [591 E]- We conclude 
that the spin-dependent part of the Id — 2p proton-neutron interaction is small. This is 
possibly due to the fact that the monopole terms ^id"/22pi/2 '^^zi^i-Vii^ weak in 
absolute value, as both their numbers of nodes and orbital momenta differ by one unit. 

From these observations, it was inferred also in Ref. [HB] that the change of p SO 
splitting {bSO{j))) between ^^Ca to '^^Ar was ascribed to a modest depletion of the Si/2 
orbit by about 0.66 protons (see also Section [s]): 

8S0{p) ^ 0.66(^-^2.,/, - y^^,2.,.J (12) 

Thus at ~ 28 the effect of the two-body SO interaction is weak. The = 21 region 
is then more propicious to study this interaction as more protons are depleted in the 
2si/2 orbit between ^''S2i and ^^Si2i. 

As for '^^Si2i, energies and spectroscopic factors of its first 7/2~, 3/2~, l/2~ 
and 5/2~ states were determined recently using the (d,p) transfer reaction in inverse 
kinematics at GANIL. They are compared to those of the isotone ^''S2i [72] in 



Figure 11 The observed p SO splitting changes by about 900 keV between ^^S and 
■^^Si while protons are removed from the 2si/2 orbit. This variation is observed in the 
major fragment of the single-particle strength, which contains correlations beyond the 
monopole terms. In order to unfold these correlations and to derive the change in 
monopole, the SM approach is used. While adjusting the monopole terms used in the 
SM calculations in order to reproduce the correlated experimental data [73j , a decrease 
of the p SO splitting from 2 MeV to 1.73 MeV is found for ^^S and an increase from 1.1 
MeV to 1.35 MeV is found for ^^Si. It follows that the change in [^3/2 — P1/2] spin-orbit 
splitting {bSO{p)) amounts to about 25%, namely (1.73 — 1.35) = 0.38 MeV compared 
to the mean value (1.73 + 1.35)/2 = 1.54 MeV. This change is mainly ascribed to a 
change in the proton 2si/2 occupancy which is calculated to be A(2si/2)=1.47 : 

bSO{p) ^ lAliY!-^^,^^^^ - V?,;^2P3/J = 380 keV (13) 

Note that by using the monopole terms (V2^"^22pi/2 ^"^^ ^25*^/22^3/2) derived in Section jsj 



and Figure ^ one finds using Equation 13 a consistent change in the p SO splitting 



of 1.47(85+170)= 370 keV between and ^^SQ The present change in SO splitting 

* In practice the monopole terms should be slightly renormalized according to the A^/'^ scaling rule 
between A ~ 48 and A ~ 36 region 
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Figure 11. Distribution of the major fragments of the single-particle strengths in the 
iV = 21 isotones. Spectroscopic factors are represented in the y axis (adapted from 
Ref. 1121). 



is exclusively due to the two-body spin-orbit force, as monopole terms involving an 
i = component, as it does here with the Si/2 orbit, does not contain any tensor 
component [30]. The only worry for this extraction of two-body term comes from the 
fact that the energy and wave function of l/2~ state in ^^S, which is bound by about 
0.5 MeV, may be influenced by the proximity of the continuum. 

Having an orbital momentum i=0, the 2si/2 orbit is located in the center of the 
nucleus. With a weak 2si/2 occupancy, ^^Si is expected to have a depleted proton central 
density [H] in the SM as well as in the MF approaches. This depletion is also seen in 
the density profile calculated using the RMF/DDME2 interaction [73] in the left part 
of Figure [121 On the other hand, the proton densities of the other = 20 isotones do 



not display any central depletion (see the middle and the right part of Figure 12). With 
this peculiarity of density depletion and isospin difference in central density (proton 
and neutron density profiles differ significantly), the ^^Si nucleus is a good candidate to 
study the SO interaction in mean-field approaches, as described below. 



6.3.3. The spin-orbit interaction in the mean- field approaches Interesting is 
to note that all theoretical approaches do not agree on the reason for the = 28 shell 
disappearance. Li et al [18j say that the ^RMF model automatically reproduces the 




Figure 12. Proton density profiles obtained with the RMF/DDME2 calculation in the 
= 20 isotones. From the left-hand picture, it is seen that the '^^Si nucleus displays 
a central depletion. By courtesy of J. P. Ebran. 

= 28 spherical shell gap because it naturally includes the spin-orbit interaction and 
the correct isospin dependence of this term'. Therefore in the RMF approach Hhere is 
no need for a tensor interaction to reproduce the quenching of the spherical N = 28 gap 
in the neutron-rich nuclei'.. The bubble nucleus ^"^Si will be used to ascertain these 
statements. 

Relativistic Mean Field (RMF) models introduced description of the nuclear spin- 
orbit interaction [76] in terms of mesonic degrees of freedom, treating the nucleons as 
Dirac particles. Using a non-relativistic reduction of the Dirac equation, the spin-orbit 
term writes: 

V;^^(r) = -[WidrPrir] + W2drPr^r'{r)]7 ■ ^ (14) 

in which Wi and W2 depend on the a, u, p meson coupling constants, and r represents 
a proton or a neutron. Note that the Wi and W2 parameters have also density 
dependence in the RMF calculations that we neglect here. Thus beyond the i ■ s 
term, the SO interaction contains a density dependence through drp{r) and an isospin 
dependence through the W1/W2 ratio. Both relativistic and non-relativistic mean-field 
(MF) approaches agree on a significant density dependence of the spin-orbit interaction. 
However, while there is moderate isospin dependence of the SO interaction in the RMF 
approach {W1/W2 — 1), non-relativistic Hartree-Fock approaches have a large isospin 
dependence {Wi/W2=2) provided by the exchange term |77j of the nuclear interaction. 
Beyond the fundamental importance to describe the spin orbit interaction, the density 
and isospin dependences of the spin-orbit interaction have important consequences to 
model nuclei close to the drip-line [78] , to predict the isotopic shifts in the Pb region jEH] 
and to describe superheavy nuclei which display a central-density depletion [80] . 

Having a significant proton-density depletion and a proton-neutron density 
asymmetry (there is no neutron central depletion), the ^^Si nucleus is a good candidate to 
study both the density and isospin dependence of the spin-orbit interaction as proposed 
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in Ref. [72]. Indeed RMF models give a large change in SO splitting by about 80% 
between ^^S and ^^Si while MF models find a more modest change of about 30% for 
A(2si/2) =1.5 [72]. 

The presently observed decrease of the SO splitting between ^''S and ^^Si by about 
25% indicates that there is a density dependence of the SO interaction, otherwise 
no change in SO splitting would have been found. This modest change is in better 
agreement with theoretical models having a large isospin dependence, such as most of 
the MF models. Though extremely preliminary, this discovery would suggest that the 
treatment of the isospin dependence of the SO interaction in the RMF approach is not 
appropriate. If true this would have important consequences to model nuclei which are 
sensitive to the isospin dependence of the SO interaction. We give here two examples. 

In some superheavy nuclei, large proton and neutron central-density depletion may 
be present. They originate from the coulomb repulsion between protons and from 
the filling of high-j neutron orbits which are located at the nuclear surface only [80] . 
Assuming no isospin dependence of the SO interaction as in the RMF approach, the 
proton and neutron central-density depletions mutually reduce the proton SO splitting 
for (£ =1, 2, and 3) orbits which are probing the interior of the nucleus [BD]. In this 
case, the Z = 114 gap, formed between the proton /5/2 and fj/2 would be significantly 
reduced to the benefit of a large shell gap at Z = 120, = 172. Conversely, a rather 
large isospin dependence of the SO interaction would not lead to strong shell gaps [HI] 
but rather to sub-shells the rigidity of which is further eroded by correlations. The 
second example deals upon nuclei close to the neutron drip-line. As they are expected 
to have a smoothened neutron density at their surface, the spin-orbit of high-j orbits is 
expected to be weakened. 

In the RMF approach the SO interaction is weaker by about up to 65% as compared 
to MF approaches when reaching the drip -line [TS]. This change between RMF and 
MF models seems to be due to the different treatment of the isospin dependence of the 
SO interaction [78]. This change in SO interaction could have important consequences 
to model the evolution of the A^ = 28 shell closure far from stability as well as other 
SO magic numbers along which the r-process nucleosynthesis occurs [82]. We note here 
that other competing effects come into play to modify shell structure when approaching 
the drip-line, as mentioned in the following Sections. 

6.4- Would other SO magic numbers 14, 50, 82 be disappearing far from stability as 
N = 28 does and if not why ? 

6. 4-1- Introduction In this Section we look whether the SO magic gaps 14, 50, 82 are 
progressively vanishing when moving far from stability in the same manner as A^ = 28 
does toward ^^Si. We start with the analogy of the evolution of the A^ = 14 shell gaps 
between ^^O and and pursue with the evolution of the SO magic number in the ^^^Sn 
nucleus at A^ = 82. We show that, surprisingly, while the same forces are involved in all 
regions, the A^ = 82 gap remains very large contrary to the A^ = 14 and A^ = 28 ones. 
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The ^®Ni nucleus, which hes in between the ^^Si and ^^^Sn nuclei, is briefly discussed as 
well. 

6.4-2. Disappearence of the = 14 and = 28 gaps far from stability The 

iV = 28 gap vanishes progressively by the removal of 6 protons from the doubly magic 
^^Ca to the deformed ^^Si nucleus. As far as tensor forces are concerned, their action 
cancels in a spin saturated valence space, i.e., when orbits with aligned and anti-aligned 
spin orientations with respect to a given angular momentum are filled. In the present 
valence space for instance, the filling of the proton fi/2 and /5/2 annihilates the tensor 
effects on the filled ^5/2 and (^3/2 orbits. As the ^^Si nucleus is a spin unsaturated system 
in protons and neutrons, tensor effects are maximized, leading to reductions of both the 
proton d3/2 — ^5/2 and neutron fi/2 — f^/2 splittings. Similar spin unsaturated nuclei 
exist in nature, such as (i) (vrpi- and ud^ filled but not vrpj^ and z/rf^), (ii) '^'^Ni (tt/i- 
and vg^ filled but not vr/^ and vgi) and (ii) ^^^Sn ij^g^ and z//i-|- filled but not Tig^ and 



vh^j. These examples are illustrated in Figure 13 We would expect that the behavior 
of all these four mirror valence nuclei is similar, but it is not the case. 

Starting with the first example, A^ = 14, it was found in Ref. [HS] that the energy 
of the first 2+ state drops by a factor of two between ^"^O (3.2 MeV [H]) and ^^C (1.588 
MeV). Then the reduction of the A^ = 14 gap, formed between the neutron ^5/2 and 
Si/2 orbits, could be ascribed to the two monopole interactions involved when removing 
2 protons from the pi/2 orbit between ^^O and ^"^C, i.e.: 

5G-(14) ^ 2(n- ,2.,/, - (15) 

We can estimate the reduction of the A^ = 14 gap by looking at similar forces in the 
valley of stability. The ordering of the 6/5/2 and Si/2 is exchanged between ^"^O (g.s. 
5/2+, 1/2+ at 870 keV) and ^^C (g.s. 1/2+, 5/2+ at 740 keV). Such a swapping of about 
1.6 MeV between the two neutron orbits means that l^ip^/jidg/J >> l^ip"/22si/2l- The 
large value of the |^ip'^/2i'i5/2 1 nionopole can be ascribed to the attractive part of the 
tensor term. Assuming that the two monopoles are of similar intensity between A = 16 
and A = 22, the A^ = 14 gap is expected to be reduced by 1.6 MeV between ^^O and 
^°C, starting from a value of about 4 MeV in ^^O (see Figure 10). This large reduction 



of the A^ = 14 spherical gap, by the removal of only two protons, can easily account for 
the decrease in 2+ energy in ^*^C. 

For the A^ = 28 shell gap, the removal of 6 protons is required to change drastically 
the nuclear structure. The effect of residual forces per nucleon is less important than 
in the lighter nuclei. It is also noticeable that the size of the A^ = 28 gap (4.8 MeV) is 
about 1 MeV larger than that of A^ = 14. These two effects, weaker monopole terms 
and increasingly larger SO splitting, could explain why the vanishing of the shell gap 
occurs further from stability at A^ = 28. 

6.4.3. Persistence of the A^ = 82 gap in ^'^^Sn On the other hand, ^^^Sn has all 
the properties of a doubly-magic nucleus with a high 2+ energy of 4.04 MeV [85], due 
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Figure 13. Schematic evolution of tlie SO shell gaps at N =14, 28, 50 and 82 as 
a function of the filling of the proton £ I orbits, with i =1, 2, 3 and 4, respectively. 
While the = 14 and iV = 28 magic gaps vanish in nuclei far from the stability valley 
(right column), the = 50 gap slightly decreases and the iV = 82 one remains very 
robust. 



to a large shell gap |3]. Large spectroscopic factors are found in the ^'^^Sn (d,p) ^^^Sn 
reaction [691 ES], which is a further confirmation that ^^^Sn is a closed-shell nucleus. 
Note that the = 82 gap is not strictly formed between the neutron hii/2 and /7/2 
orbits, as the hii/2 orbit lies below the ds/2 one for Z = 50 [87J. This arises from the 
fact that the /iii/2 — /i9/2 SO splitting is larger than the other SO splittings (79/2 — 97/2, 
fj/2 — /5/2 and o?5/2 — c?3/2 in the = 50, = 28 and = 14 regions respectively. This 
increasingly large splitting with i value is due to the T" ■ ^ term of Equation 



14 



6.4.4. The A^ = 50 shell gap in '^^Ni In the A^ = 50 case, the rigidity of the '^^Ni 
nucleus with respect to quadrupole excitations is not known so far. Based on the study 
of spectroscopic information between Z = 30 and Z = 38, the evolution of the A^ = 50 
shell gap has been studied in Ref. [M]- A rather modest reduction of about 0.55 MeV 
was proposed between Z = 38 and Z = 28. Nevertheless, as this reduction is combined 
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with that of the Z = 28 gap |65], it is not absolutely sure that Ni would remain 
spherical in its ground state. 

6.4-5. Conclusions Though similar tensor and spin-orbit forces are at play for 
N = 14, 28, 50 and 82, the corresponding shell closures do not have the same behavior 
when going far from stability. Quadrupole correlations dominate in = 14 and N = 28 
over the spherical shell gaps, while the = 82 gap remains remarkably rigid in ^'^^Sn. 
This feature likely comes from two effects : (i) the SO shell gaps become larger for higher 
i orbits and (ii) the monopole terms involved to reduce the shell gap are weakened in 
heavy nuclei because of the larger sizes of the orbits. In the next Section we explore the 
forces which come into play even further from stability. 

6.5. Which forces come into play in nuclei located below ^"^Si, "^^Ni and ^^"^Sn? 

6.5.1. Introduction Advances in radioactive beam productions offer the possibility 
to produce and study nuclei at the limit of particle stability. There, new forces are 
involved and the interaction with continuum states should progressively play a decisive 
role. We start this Section by a qualitative analysis of the nuclear forces involved at 
A^ = 28 below "^^Si. We propose that this gap is further reduced, and illustrate our 
purpose by looking at recent experimental studies which evidence the inversion between 
the /7/2 and P3/2 orbits. We generalize this mechanism to other shell gaps, A^ = 50 and 
A^ = 82, and propose astrophysical consequences for the r-process nucleosynthesis below 
i32Sn. 

6.5.2. The A^ = 28 shell gap below ^^Si It was shown in Refs. [HSl El] that the 
neutron drip line extends further A^ = 28 in the Si and Al isotopic chains, i.e., at least 
up to ''^Siso and ^^Also, respectively. If the A^ = 28 shell gap was large, a sudden drop 
in S2n would be found there and the drip-line would lie at A^ = 28 and not beyond. To 
give an example, the neutron drip line is exactly located at a magic shell in the O chain, 
where ^"^O is bound by about 3.6 MeV and ^^O is unbound by about 770 keV [BS]- So in 
this case, ^^O is a doubly-magic nucleus [90l El [92] and the spherical shell gap is large 
enough to hamper the onset of quadrupole correlations. By contrast the fact that the 
drip line extends further from stability at A^ = 28 is an indication for the presence of 
deformed nuclei which gain in binding energy due to correlations. This also hints for a 
further reduction of the A^ = 28 shell gap. This hypothesis is hard to prove in a direct 
way as it would require to determine the size of the spherical shell gap there, which is 
out of reach so far, and which would not be an observable when nuclei are deformed. 
Therefore we shall use qualitative arguments as well as an example taken from another 
region of the chart of nuclides where similar forces are present to predict the behavior 
of the N = 28 shell gap further from stability. With this in mind, we will propose some 
extrapolations to other regions of the chart of nuclides. 

Below ^^Si the change of the A^ = 28 shell gap between the I/7/2 and 2^3/2 orbits 
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Figure 14. Schematic evolution of the SO shell gaps at 28, 50 and 82, as a function of 
the filling of the proton £ f orbits, with £ =2, 3 and 4, respectively. The attractive pn 
interaction 1^5/2 — 1/7/2 is larger than the 1^5/2 — 2^3/2 one. Therefore the N = 28 gap 
is larger when the o?5/2 orbit is completely filled (left part). Conversely, when emptying 
the orbit (i.e., going very far from stability) the = 28 gap is expected to vanish 

(right part). The same mechanism is foreseen for the N =50 (82) gap when the 77/7/2 
(59/2) orbit is empty. 



is driven by the difference between the ^id?/22p3/2 ^id"/2i/7/2 i^o^^opol^ terms. On 
top of this monopole-driven effect, correlations play an important role in these nuclei. 
For these two monopole terms, the proton and neutron spin orientations are aligned 
with the orbital momentum, and their difference in angular momentum is one unit of 
h (then the main part of the proton-neutron interaction comes from the central term 
of the nuclear force). On the other hand, the numbers of nodes in their wave functions 
are different. The fact that the I/7/2 and 1^5/2 wave functions have the same number of 
nodes leads to a larger radial overlap of the wave functions, and a larger monopole term 
^id"/2i/7/2 compared to ^1^^/22^3/2 • Consequently when the 1(^5/2 orbit is completely 
filled, the I/7/2 orbit has a larger binding energy than that of the 2^3/2 one. Conversely, 
when removing the six protons from the 1^5/2 orbit, the neutron = 28 gap further 
shrinks. Note that the calculations of Ref. [32] predict a decrease by 1.6 MeV down to 
■^^028- Of course this nuclei lies beyond the drip-line, but this gives the monopole trend 
for the evolution of the N = 28 gap below ^^Si (see Figure 14). 

In addition to this monopole effect, the nuclear-potential well will become more 
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diffuse when moving towards the drip-hne. It was proposed by Dobaczweski et al [93] 
that this would bring the low i orbits more bound as compared to the high £ ones, which 
could be viewed schematically as a reduced £^ term in the Nilsson potential. Then, an 
inversion of the two negative-parity orbits (/ and p) could also happen, as suggested by 
several experimental works described below. 

While the z/l/7/2 orbit lies below the v2p^i2 orbit in most of the nuclei in the chart 
of nuclides (giving rise to the = 28 shell gap), their ordering is reversed as soon as 
the 7r(i5/2 shell starts to empty, i.e., for Z < 14. The feature is well documented in 
the = 15 and = 17 isotones, not far from the stability valley. Indeed, the first 
negative-parity state of i4^^Sii5_i7 has P = 7/2^ while the 3/2~ state is located above. 
On the other hand, the situation is reversed in i2Mg and loNe isotopes, as described 
now. In the ^^Neis nucleus [Mj, the 2^3/2 orbit lies about 1 MeV below the I/7/2 orbit. 
This was determined by using the ^'^Ne(d,p)^^Ne transfer reaction in inverse kinematics 
at the GANIL facility. Similar results were recently obtained at GANIL in ^''Neiy, its 
first excited state at 765 keV having I'^ = 3/2~ [95], the 7/2~ state being about 1 MeV 
above. Note that the assignment of the 3/2~ state was first proposed in Refs. [961 197] . 
In the neighboring isotone ^^Mg, the second excited state at 1095 keV was observed from 
the study of the /3 decay of ^^Na. It could not be explained in the framework of shell- 
model calculations using the sd space, and consequently it was assigned as a /'^ = 3/2~ 
intruder state [98]. To our knowledge, two other experimental works are suggesting 
this inversion between the 1/ and 2p orbits. By using Coulomb break-up technique at 
RIKEN, Nakamura et al [99j found that the ground-state wave function of ^^Ne was 
dominated by an £ = 1 component. Minomo et al |100j propose a contribution of a P3/2 
neutron using antisymmetrized molecular dynamics calculations. Moreover Wimmer et 
al argued that excitations to the 2^3/2 (rather than the I/7/2) are required to account 
for the configuration of the Of 2 states in ^^Mg determined through two-neutron transfer 
reaction at CERN/ISOLDE [lOlj . In these examples, in which the proton 1^5/2 orbit is 
not yet fully filled, the = 28 gap does not exist and the 1^2^3/2 orbit even lies below 
the plfT/2 orbits. Conversely, this gives a further indication that the proton 1^5/2 has 
to be filled to bind the 1^1/7/2 orbit enough to create a A^ = 28 gap which is present 
in heavier nuclei. Note that, from a theoretical point of view, this inversion was also 
proposed by Utsuno et al in 1999 |102j . 



6.5.3. Extrapolations to other shells A^ = 50 and A^ = 82 We can reasonably 
assume that the present mechanism for the reduction of the A^ = 28 gap is robust, 
as it comes from the radial part of the wave functions and from the effect of a more 
diffused potential well at the drip-line. As far as the first effect is concerned, it is also 
observed from the properties of the bare proton-neutron forces |103j . We therefore also 
expect the reduction of the A^ = 50 gap below ^^Ni and the reduction of the A^ = 82 
gap below ^^^Sn, as protons are removed from the I/7/2 and 1(79/2 orbits, respectively 
(see Figure 14). While the A^ = 82 gap is still present in the ^^^Sn nucleus, it will 
progressively be reduced further from stability as ten protons are removed from the 
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1(79/2 orbit to reach ^^^Zr. The onset of quadrupole correlations across the = 82 shell 
gap would occur if the neutron gap is reduced enough. This would bring a smoothening 
in the S2n trend, as observed in Figure [2] for the N = 28 isotones around ''^S. This has 
potentially important consequences on the location of the (n, 7) — (7, n) equilibrium in 
the rapid neutron-capture process (r-process). This equilibrium occurs when neutron 
captures compete with photo-disintegrations. At this point, further neutron captures 
cannot occur as nuclei are immediately destroyed. Therefore the explosive process is 
stalled for a moment at such waiting-point until /3-decay occurs. When the drop in 

value is abrupt (for a large shell gap), the waiting-point nuclei are found at the 
closed shell. Then these nuclei are the main genitors of stable r-elements when decaying 
back to stability. When the shell gap is reduced, the location of r-progenitors is more 
extended and shifted to lower masses in a given Z chain, changing the fit of the r-process 
abundance curve accordingly (see for instance Ref. [82]). 

To summarize this Section, we propose that both from the nuclear two-body forces 
and from the mean-field point of views a further reduction of the = 28 shell gap will 
be occurring below '^^Si. By analogy, a reduction of the = 50 and = 82 shell gaps 
is anticipated, with possible important consequences for the r-process nucleosynthesis 
in the later case. 

6.6. Are proton-neutron forces significantly modified when approaching the drip line ? 

6.6.1. Introduction In the previous Sections we assumed that the two-body proton- 
neutron forces (or the monopole terms) remain of similar intensity when reaching the 
drip line. The SM description does not allow a self-consistent change of the monopole 
terms as a function of binding energy. However the assumption of fixed monopole 
terms in a wide valence space is at best crude and probably wrong. Indeed at the 
neutron drip-line, valence protons are very deeply bound and their wave functions are 
well confined. On the other hand, neutrons, which are close to be unbound, have more 
dilute wave functions and interactions with continuum states occur |104] . It follows that 
the radial part of the wave functions of the valence proton and neutron have weaker 
overlap, leading to a reduced effective interaction. One could question by how much 
this interaction is reduced ? 

The 9F isotopic chain is a target choice to study changes in proton-neutron forces 
when approaching the drip line. Indeed the drip-line occurs at A^ = 16 in the sO chain, 
i.e., at ^^O, while the ^^F22 nucleus is still bound. Thus the addition of a single proton 
allows to bind up to 6 neutrons more. It follows that, though occurring between nuclei 
showing large asymmetry in neutron and proton binding energy, the corresponding 
proton-neutron forces should be large enough to bind 31F2II Up to 29F20, the same 
proton- neutron d^j2 — d^/2 is expected to intervene. This strong monopole force is well 
known to produce significant changes closer to stability in the shell evolution of the 

Beyond proton-neutron interactions, we remind here the important role of neutron-neutron forces as 
weh. 
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= 16 and = 20 gaps. How is this force modified far from stability ? The recent 
determinations of the binding energies of the (unbound) ^^Fig |105j and the bound 
^^F2o |106j offer the possibihty to start answering these questions. However error bars 
on binding energies are still a bit too large at present to draw firm conclusions. 

In this Section we propose to study the spectroscopy of the l^Fn to study the 
proton-neutron ^5/2 — 1^3/2 interaction when approaching the drip line. We first explain 
the motivation of studying this nucleus and subsequently use the new spectroscopic 
information required to study the proton-neutron force. 



6. 6.2. Experimental study of g^F The 9^Fi7 nucleus is a suitable choice for studying 
the proton-neutron 0^5/2 — interaction for at least three arguments: (i) it is bound 
by only 0.80(12) MeV, (ii) as the first excited state in ^^Q lies at 4.47 MeV [9l], ^^F can 
be viewed as a closed ^"^O core plus a deeply hound proton in the d^/2 orbit (S'p(^^F)~ 
-15 MeV) and an unbound neutron in the ^3/2 orbit (S'„(2^0)~ -F0.77 MeV (SH]), (iii) 
the TTd^/2 and i^c?3/2 orbits are rather well separated from other orbits which limit 
correlations of pairing and quadrupole origin. Binding energies BE(^^F)j of the full 
multiplet of J = 1 — 4+ states arising from the 71^5/2 ® i'd-i/2 coupling in '^'Pn are 
needed to determine the role of the coupling to the continuum on the mean 'Kd^j2i'd^/2 
interaction energy (mt), defined in Equation (16), as well as on the residual interaction 
which lifts the degeneracy between the components of the multiplet int{J) defined in 
Equation (17). The mean interaction energy (int) writes: 

^^^^^^i2J^_ll^ (16) 

In this Equation int{J) term expresses the difference between the experimental binding 
energy of a state J in ^^F {BE(^^F)j) and that of the ^^Q+lp+ln system, BE(^^Ffree), 
in which the valence proton and neutron do not interact. It writes: 

mt(J) = BEC^F)j - BEC^Ffree), (17) 

where 

BE^^Ffree) = BE{^^0)o+ + BE{nd5/2) + BE{iydy2) (18) 
Assuming that 

BE{nd,/2) = BEC'>F),/2+ - BEC'0)o+ (19) 

and 

BE{vd^/2) = 5^(2^0)3/2+ - BEC''0)o+ (20) 

we obtain: 

BE{^'Ffree) = BEC'F),/2^ + BE{^'>0),/2^ - BE(^'0\^ (21) 

Following the particle-particle coupling rule, the values of int (J) should display a 
parabola as a function of J in which |m)f:(l)| and |mt(4)| are the largest, as proton 
and neutron maximize their wave-function overlap for these two values of total angular 
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Figure 15. Schematic view of the low-energy configurations in ^^F, which could be 
approximated as a core of on top of which a deeply bound proton ^5/2 and an 
unbound neutron are added. This could bring, after having treated the role of 
correlations, information on the proton-neutron — ^3/2 interaction close to the 
drip-line. In the bottom part, values of int(J) extracted from the SM calculations are 
compared to the presently known experimental values for J'^ = 1"*" — 4"*". 



momentum, leading to the largest attractive interactions. The values of int{J) obtained 
from shell model calculations (see bottom part of Figure 15) indeed form a parabola, 
the \int{3)\ value being the lowest. 

The atomic masses of 26,25p 24q t^qj-q determined by Jurado et al They were 
used to determine the binding energies BE values used in the previous Equations. The 
5^(^^0)3/2+ - BE(^'^0)o+=+0.77 MeV difference is derived from the work of Hoffman 
et al [89] . The spin of the ground state of ^^F was determined to be f = 1"*" by Reed 
et al |1U7] from the fact that f = 0,2+ states were populated in its (3 decay. The 
first-excited state of ^^F has been recently identified at 657(7) keV excitation energy by 
Stanoiu et al |108j at GANIL using in-beam 7-ray detection combined to the double 
fragmentation method. Moreover, an unbound state lying 270 keV above the neutron 
emission threshold has been proposed by Franck et al following a nucleon-exchange 
reaction at NSCL |109] . This state is likely to be the I'^ = state. Finally a long-lived 
J = 4+ isomer, T1/2 = 2.50(5)ms, has been discovered at 642 keV by Lepailleur et al 
at GANIL [llOj . It decays by an M3 internal transition to the J = ground state 
(86%), and by beta-decay to ^^Ne or beta-delayed neutron emission to ^^Ne. Gathering 
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the measured binding energies of the J = 1 — 4"'" muhiplet in g^F in Figure 
that the components of the muhiplet are experimentaUy more compressed than the SM 
predictions, whatever the effective interaction, USDA or USDB . This would be due 
to a weaker residual interaction 71^5/2 z^c?3/2 close to the drip line. It would be interesting 
to compare these findings to calculations taking into account the effect of continuum. 
Note that the present conclusions rely strongly on atomic mass determination of the g^F 
ground state, as well as the resonance at ~270 keV we have assigned to be the I'^ = 3"*" 
state of the multiplet. These two results deserve to be confirmed. 

7. General conclusions 

The evolution of the = 28 shell closure has been investigated far from stability, 
showing a reduction of the spherical shell gap. It was found that a progressive onset of 
deformation occurs after the removal of several protons. Starting from a rigid-spherical 
2oCa nucleus, fgAr has a vibrational character, f^S exhibits a shape coexistence, ^481 is 
likely to be oblate and f^^g is expected to be prolate. Adding the triaxiality in the 
^^Ar nucleus, a wealth of nuclear structures can be explored in this mass region. 

This structural evolution has been probed using various experimental techniques, 
such as in-beam 7-ray spectroscopy, transfer. Coulomb excitation, g factor measurement; 
to quote a few. The evolution of the sd proton and fp neutron orbits was probed, 
from which important conclusions were derived. Starting with proton orbits, a total 
disruption of the Z = IQ sub-shell, formed between the Si/2 and ^3/2 orbits, occurs from 
A^ = 20 to A^ = 28. In parallel a reduction of the Z = 14 gap (between the and Si/2 
orbits) seems to be present as well. This global shrink between proton states favors the 
development of pairing and quadrupole excitations at A^ = 28 as compared to A^ = 20, 
where the magicity persists down to ^481. As for neutrons orbits, reductions of the / and 
p SO splittings as well as the A^ = 28 shell gap (by about 25%) were proposed between 
"^"^Ca and ^^Si. This reduction is large enough to trigger the development of quadrupole 
excitations across the N = 28 gap. In all, both the proton and neutron shell gaps are 
reduced in the spin-unsaturated ^^Si nucleus, mutually enhancing the collectivity there. 

These studies have in addition shown the many facets of nuclear force that 
can be explored, such as the central, spin-orbit, tensor and three-body terms. It 
also triggered the question of proton- neutron forces below ^^Si and more generally 
when approaching the drip-line. These points were addressed in the last part of the 
manuscript. Conclusions are as follows: 

• Due to three-body effects, the N = 28 gap grows from A^ = 20 to A^ = 28 by 
about 2.7 MeV. As similar forces seem to be at play in the other SO shell gaps, an 
empirical rule was proposed to predict the increase of the A^ = 14, A^ = 50, A^ = 82 
shell gaps when the i t orbits (with i =2,3,4) are completely filled. Following this 

ft The atomic mass of ^^F measmed in [4] was shifted to a larger binding energy (downward) to take 
into account the probable contamination of 40% in isomeric state. 



15|[tt we find 
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rule, a new = 90 sub-shell gap should be present in ^ Sn, leading to interesting 
consequences for the r-process nucleosynthesis. 

• The roles of the central, tensor and spin-orbit components of the nuclear force were 
examined in the framework of the shell model using various experimental data. 
The bubble nucleus ^''Si has been used to probe the reliability of the density and 
isospin dependence of the SO interaction in mean-field models. It was derived that 
the spin-orbit interaction is indeed density dependent, but also seems to have an 
important isospin dependence that is not present in RMF approaches. This feature 
should have significant consequences to model the structure of superheavy and 
neutron-rich nuclei, which exhibit central density depletion and surface diffuseness, 
respectively. 

• Below ^^Si, central forces further reduce the = 28 gap. The same situation is 
expected to occur for the = 50 and = 82 shell gap below ^'^Ni and ^^^Sn, 
respectively. The reduction of the A^ = 82 gap would have extremely important 
consequences for the r-process nucleosynthesis. 

• When reaching the neutron drip-line, proton-neutron forces may be weakened by the 
fact that neutrons are no longer confined in the nuclear potential well. The study 
of ^^F was proposed to test such effects. It was found that both the monopole 
and the residual proton-neutron interactions are weakened close to the drip-line 
as compared to the values used for nuclei close to stability. This very preliminary 
study is foreseen to be a benchmarking case for shell-model approaches or coupled 
cluster models which treat the influence of the continuum. 

More generally, in the present work we have explored in a semi-quantitative manner 
several facets of the nuclear force which lead to modification of shell structure far from 
stability. Most of these properties were derived from two-body effective interactions 
(including monopoles) in the shell-model framework. The effective interactions are built 
from realistic bare interactions to which medium effects are added. These interactions 
are subsequently fitted to experimental data to include effects from missing terms such 
as the three-body components. The search for which parts of the nuclear interaction give 
rise to significant shell modifications was started a long time ago (e.g. in |1121I113[I114] ). 
However very recent works manage to point out which properties of the bare nuclear 
force are preserved in the nucleus and which are modified. To quote a few it was found 
in Refs. [311 EO] that the tensor force is almost kept constant from the bare forces to the 
nucleus, while the central part of the interaction varies a lot. Moreover some theoretical 
works can now model nuclear systems with few nucleons or close to a core nucleus 
starting from realistic interactions or/and treat the interaction with the continuum. 

As for experiments, important results are to come soon in the A^ = 28 region, 
such as the B(E2) value and the identification of higher- lying states in ^^Si, as well as 
the spectroscopy of neutron- rich Mg isotopes when approaching A^ = 28. Tremendous 
progresses have been achieved over the last 20 years since the first suggestion that the 
"^^S nucleus was deformed, based on its very short lifetime and small neutron delayed 
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emission probability |115] . to the generalization of nuclear properties in all SO magic 
nuclei. All this was possible with the increase of experimental capabilities as well as 
regular exchanges between experimentalists and theoreticians. 
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